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We consider direct diphoton production in hadron collisions. We compute the next-to-next-to- 
leading order (NNLO) QCD radiative corrections at the fully-differential level. Our calculation 
uses the qr subtraction formalism and it is implemented in a parton level Monte Carlo pro- 
i <""| gram, which allows the user to apply arbitrary kinematical cuts on the final-state photons 

O-i an d the associated jet activity, and to compute the corresponding distributions in the form 

of bin histograms. We present selected numerical results related to Higgs boson searches at 
the LHC, and we show how the NNLO corrections to diphoton production are relevant to 
r*j . understand the main background of the decay channel H — > 77. 



J> 1 1 Introduction 

(NT 

Diphoton production is a relevant process in hadron collider physics. It is both a classical 
^r) " signal within the Standard Model (SM) and an important background for Higgs boson and new- 

physics searches. Recent results from the LHC indicates that the Higgs boson mass mjj must 
be low (114 GeV< mn < 130 GeV), and thus the preferred search mode involves Higgs boson 
production via gluon fusion followed by the rare decay into a pair of photons. We are interested 
in the process pp — > 77X, which, at the lowest order in perturbative QCD, occurs via the quark 
annihilation subprocess qq —> 77. The QCD corrections at the next-to-leading order (NLO) in 
^ the strong coupling as involve the quark annihilation channel and a new partonic channel, via 

the subprocess qg — > 779. These corrections have been computed and implemented in the fully- 
differential Monte Carlo codes DIPHOX,^ 2gammaMcl21 and MCFmP A calculation that includes the 
effects of transverse-momentum resummation is implemented in RESBOS.^ At the next-to-next-to- 
leading order (NNLO), the gg channel starts to contribute, and the large gluon-gluon luminosity 
makes this channel sizeable. Part of the contribution from this channel, the so called 60a; 
contribution, was computed long ago^ and its size turns out to be comparable to the lowest-order 
result. Besides their direct production from the hard subprocess, photons can also arise from 
fragmentation subprocesses of QCD partons. The computation of fragmentation subprocesses 
requires (poorly known) non-perturbative information, in the form of parton fragmentation 
functions of the photon. The complete NLO single- and double-fragmentation contributions are 
implemented in DIPHOX^The effect of the fragmentation contributions is sizeably reduced by the 
photon isolation criteria that are necessarily applied in hadron collider experiments to suppress 
the very large irreducible background (e.g., photons that are faked by jets or produced by hadron 
decays). The standard cone isolation and the 'smooth' cone isolation proposed by Frixione^ 
are two of these criteria. The standard cone isolation is easily implemented in experiments, 
but it only suppresses a fraction of the fragmentation contribution. The smooth cone isolation 
(formally) eliminates the entire fragmentation contribution, but its experimental implementation 



is still under study^ However, it is important to anticipate (work to appear), that in some 
kinematical regions (e.g for Higgs boson searches), the standard cone and the Prixione isolation 
criteria give basically the same theoretical answer § 



2 Diphoton production at NNLO 

We consider the inclusive hard-scattering reaction hi + hi — > 77 + X, where the collision of 
the two hadrons, h\ and hi, produces the diphoton system F = 77 with high invariant mass 
M 77 . The evaluation of the NNLO corrections to the this process requires the knowledge of 
the corresponding partonic scattering amplitudes with X = 2 partons (at the tree level , ^j) 
X = 1 parton (up to the one-loop leveJ^J) and no additional parton (up to the two-loop level^J) 
in the final state. The implementation of the separate scattering amplitudes in a complete NNLO 
(numerical) calculation is severely complicated by the presence of infrared (IR) divergences that 
occur at intermediate stages. The subtraction formalism^] j s a method that handles and 
cancels these unphysical IR divergences up to the NNLO. The formalism applies to generic 
hadron collision processes that involve hard-scattering production of a colourless high-mass 
system F. Within that frameworkp^ the corresponding cross section is written as: 



da (N)NLO — ^-(N)NLO ® d^LO + 



1 F+jets , CT 
aa {N)LO aa (N)LO 



(1) 



where da,^}^ represents the cross section for the production of the system F plus jets at 

(N)LO accuracy §, and da9£ LO is a (IR subtraction) counterterm whose explicit expression^ 
is obtained from the resummation program of the logarithmically-enhanced contributions to qr 
distributions. The 'coefficient' T~L F N ^ NLO , which also compensates for the subtraction of da9^ LO , 
corresponds to the (N)NLO truncation of the process-dependent perturbative function 

7T \ 7T / 

The NLO calculation of da F requires the knowledge of T~L F ^\ and the NNLO calculation also 
requires T-L F ^ 2 \ The ge neral s tructure of T~L F ^ is explicitly known^^ exploiting the explicit 
results of T~L F ^ for Higgs^^^ and vector boson^l production we have generalized the process- 
independent relation of Ref?^ to the calculation of the NNLO coefficient % F(2) . 

3 Quantitative results 

We have performed our fully-differential NNLO calculation^^ of diphoton production according 
to Eq. dU). The NNLO computation is encoded in a parton level Monte Carlo program, in 
which we can implement arbitrary IR safe cuts on the final-state photons and the associated 
jet activity. We concentrate on the direct production of diphotons, and we rely on the smooth 
cone isolation criterionP Considering a cone of radius r = yj (At/) 2 + (A0) 2 around each photon, 
we require the total amount of hadronic (partonic) transverse energy Et inside the cone to be 
smaller than ET ma x(i"), 

(1 — COS T \ ^ 
1 T> ' ( 3 ) 
1 — cos R J 

where p7p is the photon transverse momentum; the isolation criterion Et < Ermax(, r ) has to 
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be fulfilled for all cones with r < R. We use the MSTW 2008 10 sets of parton distributions, 



"The use of the same parameters in both criteria is understood. 

b In the case of diphoton production, the NLO calculation of da"^^ ts was performed in ReJJ^ 



with densities and as evaluated at each corresponding order, and we consider Nf = 5 mass- 
less quarks/antiquarks and gluons in the initial state. The default renormalization (/xr) and 
factorization (hf) scales are set to the value of the invariant mass of the diphoton system, 
fJ-R = = M 77 . The QED coupling constant a is fixed to a = 1/137. 
To present some quantitative results, we consider diphoton production at the LHC (y/s = 




za 140 160 160 "-° 
M„ (GeV) A0. 



Figure 1: Left: Invariant mass distribution of the photon pair at the LHC (yfs = 14 TeV): LO (dots), NLO 
(dashes) and NNLO (solid) results. We also present the results of the box and NLO+box contributions. The inset 
plot shows the corresponding ¥L-factors. Right: Diphoton cross section as a function of the azimuthal separation 
of the two photons. Data from CMtP^ (yfs = 7 TeV) are compared to the NNLO calculatioi^r^ 



14 TeV). We apply typical kinematical cuts used by ATLAS and CMS Collaborations in their 
Higgs search studies. We require the harder and the softer photon to have transverse momenta 
^harder y Q e y anc [ p^iter -> 25 GeV, respectively. The rapidity of both photons is restricted 
to \y~\ < 2.5, and the invariant mass of the diphoton system is constrained to lie in the range 
20 GeV < M 77 < 250 GeV. The isolation parameters are set to the values e 7 = 0.5, n = 1 
and R = 0.4. We observe^l that the value of the cross section remarkably increases with the 
perturbative order of the calculation. This increase is mostly due to the use of very asymmetric 
(unbalanced) cuts on the photon transverse momenta. At the LO, kinematics implies that the 
two photons are produced with equal transverse momentum and, thus, both photons should 
have pj, > 40 GeV. At higher orders, the final-state radiation of additional partons opens a new 
region of the phase space, where 40 GeV > p?f > 25 GeV. Since photons can copiously be 

i-i yi 

produced with small transverse momentum 1 ^ the cross section receives a sizeable contribution 
from the enlarged phase space region. This effect is further enhanced by the opening of a new 
large- luminosity partonic channel at each subsequent perturbative order. In Fig. Q] (left) we 
compare the LO, NLO and NNLO invariant mass distributions at the default scales. The inset 
plot shows the K-factors defined as the ratio of the cross sections at two subsequent perturbative 
orders. We note that ~K NNLO / NLO [ s sensibly smaller than K NLO / LO , and this fact indicates an 
improvement in the convergence of the perturbative expansion. We find that about 30% of the 
NNLO corrections is due to the gg channel (the box contribution is responsible for more than 
half of it), while almost 60% still arises from the next-order corrections to the qg channel. 

Recent results from the LHC H2J2Q] an( ^ ^ e Xevatron^B show some discrepancies between the 
data and NLO theoretical calculations of diphoton production. Basically, discrepancies were 
found in kinematical regions where the NLO calculation is effectively a LO theoretical descrip- 
tion of the process. Such phase space regions! are accesible at NLO for the first time, due to the 



c Away from the back-to-back configuration. 



final-state radiation of the additional parton@ Figured] (right) shows a measurement by CMSP 
of the diphoton cross section as a function of the azimuthal angle A0 77 between the photons. 
The data are compared with our NLO and NNLO calculations?^ The acceptance criteria used 
in this analysis (Vs = 7 TeV) require of: p^ avdeT > 23 GeV and ^ ofter > 20 GeV. The rapidity 
of both photons is restricted to |y 7 | < 2.5, and the invariant mass of the diphoton system is 
constrained to be M 77 > 80 GeV. The isolation parameters are set to the values e 7 = 0.05, n = 1 
and R = 0.4. We note that the CMS data are selected by using the standard cone isolation 
criterion and the constraint in Eq. (|3j) is applied only to the cone of radius r = R. Since the 
smooth isolation criterion used in our calculation (we apply Eq. ([3]) for all cones with r < R) is 
stronger than the photon isolation used by CMS, we remark that our NLO and NNLO results 
cannot overestimate the corresponding theoretical results for the CMS isolation criterion. The 
histograms in Fig. Q] (right) show that the NNLO QCD results remarkably improve the theoret- 
ical description of the CMS data throughout the entire range of A^ 77 . 

The results illustrated in this contribution show that the NNLO description of diphoton produc- 
tion is essential to understand the phenomenology associated to this process, and therefore, the 
NNLO calculation is a relevant tool to describe the main background for Higgs boson searches. 
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